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ABSTRACT 

We present the first resolved images of the eclipsing binary (3 Lyrae, obtained 
with the CHARA Array interferometer and the MIRC combiner in the H band. 
The images clearly show the mass donor and the thick disk surrounding the mass 
gainer at all six epochs of observation. The donor is brighter and generally ap- 
pears elongated in the images, the first direct detection of photospheric tidal 
distortion due to Roche-lobe filling. We also confirm expectations that the disk 
component is more elongated than the donor and is relatively fainter at this wave- 
length. Image analysis and model fitting for each epoch were used for calculating 
the first astrometric orbital solution for (3 Lyrae, yielding precise values for the 
orbital inclination and position angle. The derived semi-major axis also allows 
us to estimate the distance of (3 Lyrae; however, systematic differences between 
the models and the images limit the accuracy of our distance estimate to about 
15%. To address these issues, we will need a more physical, self-consistent model 
to account for all epochs as well as the multi-wavelength information from the 
eclipsing light curves. 
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Introduction 



Interacting binaries are unique testbeds for many important astrophysical processes, 
such as mass and momentum transfer, accretion, tidal interaction, etc. These processes 
provide information on the evolution and properties of many types of objects, including 
low-mass black holes and neutron stars (in low-mass X-ray binaries), symbiotic binaries, 
cataclysmic variables, novae, etc. Although these types of objects are widely studied by 
indirect methods such as spectroscopy, radial velocity, and sometimes eclipse mapping, very 
few of them have been directly resolved because they are very close to each other and far 
away from us. Thus, directly imaging interacting binaries, although very challenging, will 
greatly help us to improve our understanding of these objects. 

The star (3 Lyrae (Sheliak, HD 174638, HR 7106, V = 3.52, F=3.35) is a well known 
i nteracting and eclipsing binary that has been wi dely studied sinc e its discovery in 1784 
( )Goodrickdll785l ) . According to the current picture (lHarmanedl2002l ) . the system consists of 
a B6-8 II Roche-lobe filling mass-losing star, which is generally denoted as the donor or the 
primary, and an early B type mass-gaining star that is generally denoted as the gainer or the 
secondary. The donor, which was initially more massive than the gainer, has a current mass 
of about 3 M , while the gainer has a mass of about 13 M . It is thought that the gainer is 
completely embedded in a thick accretion disk with bipolar jet-like structures perpendicular 
to the disk, which creates a light-scattering halo above its poles. The orbit of the system 
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is highly circular (lHarmanec fc Scholzl Il993l ) , and is very close to edge-o n (ILinnell 
Recent RV study on the ephemeris of the system gives a period of 12. d 94 (lAk et al. 
The period is increasing at a rate of ~ 19 sec per year due to the high mass transfer rate, 
2 x 1O" 5 M yr^ 1 , of the system. 

The primary eclipse of the light curve (i.e., at phase 0) corresponds to the eclipse of 
the donor. In the UBV bands, the surface of the donor is brighter than that of the gainer, 
and therefore the primary minimum is deeper than the seconda ry mi nimum. At longer 
wavelengths, however, the studies of iJameson fc Longmord (119761 ) and IZeilik et al.l (119821 ) 
suggest that the relative depth of the secondary minimum in the light curve gradually deepens 
and becomes deeper than the primary minimum at wavelengths longer than 3.6//m. 

Light curve stu dies and theoretical models have shown that, at the distance of 296pc 
(Ivan Leeuwenl 120071 ). the estimated separation of the binary is only 0.92 milli-arcsecond 
(hereafter mas, 58.5 Rq). The angular diameter of the don or is ^0.46 mas (29.4 R ), an d 
the disk surrounding the gainer is only ~ 1 mas across (e.g., lLinnellll2000l ; lHarmaned 120021 ) . 



The goal of directly imaging (3 Lyr, therefore, requires the angular re s olutio n only achievable 
by today's long-baseline interferometers. Recently, ISchmitt et al.l (120081 ) used the NPOI 
interferometer to image successfully the Ha emission of (3 Lyr, an update to the pioneering 
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work of lHarmanec et al.l (119961 ) . Also, radio work using MERLIN found a ne bula surrounding 
the secondary but could not resolve its bipolar shape (jUmana et al.ll2000l ). Despite recent 
progress, the individual objects of the system have not been resolved yet, putting even a 
simple astrometric orbit beyond our reach. 

In this study, we present the first resolved images of the j3 Lyr system at multiple phases, 
obtained with the CHARA Array and the MIRC combiner. We give a brief introduction to 
our observations and data reduction in §2J We present our aperture synthesis images with 
simple models in §31 In §5]we discuss our astrometric orbit of (3 Lyr and we give the outlook 
for future work in Sj5l 



2. Observations and data reduction 



Our observations were conducted at the Georgia State University (GSU) Center for 
High Angular Resolution Astronomy (CHARA) interferometer array along with the MIRC 
instrument. The CHARA array, located on Mt. Wilson and consist ing of six 1-meter tele- 



scope s, is the longest optical/IR interferometer array in the world (jten Brummelaar et al. 



20051 ) . It has 15 baselines ranging from 34m to 331m, providing resolutions up to 
at H band and ~ 0.7 mas at K band. 



^0.5 mas 



The Michigan Infra-Red Combiner (MIRC) was used here to combine 4 CHARA tele 
scopes together for true interferometric imaging in H band, providing 6 visibi lities, 4 closure 



phases and 4 triple amplitudes simultaneously in 8 narrow spectral bands (see iMonnier et al. 



20041 . 120061 . for details). Specifically, the (3 Lyr system was observed on 6 nights in 2006 and 



2007 using an array configuration optimized for good ima ging (equal Fourier co verage in all 
directions) and following standard observing procedures (IMonnier et al.l 120071 . M. Zhao et 
al. 2008, in preparation). A typical baseline coverage of our observations is shown in Figure 
[TJ In short, we observed our target along with two or three calibrators on each night and a 
complete observing log is listed in Table [TJ 



The data reduction process follows the pipeline outlined by IMonnier et al.l (120071 ). In 
brief, after frame co-adding, background subtraction and Fourier transform of the raw data, 
fringe amplitudes and phases are used to form squared-visibilities and triple products. Pho- 
tometric calibrations are estimated using shutter matrix measurements and partial beam 
chopping. Finally, calibrators with known sizes (see Table [T]) are used to calibrate the drifts 
in overall system response before obtaining final calibrated squared-visibilites and complex 
triple amplitudes. 
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3. Synthesis imaging and modeling 

For imaging with optical interferometry data, we employed two i ndependent applica- 



tions: "Markov-Chain Imager for Optical I nterferometry (MACIM)" (jlreland et al. 



and the Maximum- Entropy-based BSMEM (IBuscherl Il994l ). Further description and a de- 



20061 ) 



tailed comparison of these algorithms on simulated data appear in lLawson et al.l (12004 



20061 ). Both of these algorithms benefit from use of prior information, generally based on 
lower resolution data. For {3 Lyr we began each image reconstruction with a two-component 
Gaussian model which mainly acts to limit the field-of-view of the image. The final images 
do not resemble the priors except in general extent, i.e., the final positions, relative sizes, 
and relative brightnesses are not dependent on the priors. The final reconstructed images 
from both methods are shown in Figure [2j The MACIM and BSMEM images are consistent 
with each other, although they use very different algorithms, giving confidence to the image 
fidelity. Any differences, such as the more Gaussian shapes for BSMEM compared to the 
more "flat-top" profiles for MACIM, illustrate the limitations of our dataset. We present 
here results from both algorithms in lieu of image "error bars" which are notoriously difficult 
to define in aperture synthesis imaging. 

The six epochs span all phases of the orbit, changing from middle primary eclipse 
(phase=0.035) to nearly maximum elongation (phase=0.210 and 0.828), and secondary 
eclipse (phase=0.438). The system is well resolved into two separate components at phases 
close to the maximum elongation. Since the primary eclipse is still the deeper one in the 
if -band (Jameson & Longmore 1976), we can conclude that the object with higher surface 
brightness is the mass donor star (i.e., the component moving from left to right in the 2007 
July sequence). The donor is partially resolved and appears elongated at all epochs except at 
phase 0.035 when it is blocked by the disk, directly confirming its Roche-lobe filling picture. 
The thick disk surrounding the gainer is also resolved and appears elongated. At the first 
epoch (phase = 0.035), we see mostly the emission from the disk superposed with a small 
amount of light from the poles of the donor. 

We can extract further information by constructing a simple two-component model to 
determine the separation and position angle for each epoch. We assume the donor and gainer 
can be modeled as uniform ellipses. Other models, such as two truncated gaussian ellipses, 
a rain- drop- shaped Roche- lobe filling star with a truncated gaussian disk, etc., were also 
considered and gave equivalent results due to limited resolution. Therefore, for simplicity 
and to minimize the degrees of freedom of the model, uniform ellipses are adopted. The free 
parameters in the models are: the semi-major and semi-minor axes of the two components, 
their individual position angles, the separation and position an gle of the system , and the flux 



ratio of the donor and the disk. We used ephemeris data from lAk et al.l (120071 ) to fix which 
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component was in front during modeling. Due to degeneracies in the separation and the 
dimensions of the blocked component when the two are overlapping with each other, the size 
of the blocked component is fixed to the average from the two separated epochs, 2007Jul04 
and 2007Jull2. The best-fit models for all epochs are presented in the third column of Figure 
|2]and the resulting positions and total flux ratios from the models are listed in Table El along 
with the results obtained from the image centroids where separating the two components is 
possible. Errors of the positions are estimated from the x 2 surfaces of each parameter where 
Ax 2 = 1, or from the scatter in fits within each night, whichever is bigger. 

The models confirm that the smaller and more circular component, i.e., the donor, has 
higher surface brightness and total flux than the more elongated disk around the gainer. The 
ellipse size of the donor from the models, when averaged over all the epochs, is 0.62 ± 0.16 
mas along the major axis, and 0.52 ± 0.14 mas along the minor axis, which confirmed the 
i mages that the d onor is elongated, but slightly larger than that from the theoretical models 
jHarmanec!l2002h . The averaged size of the disk surrounding the gainer is 1.04 ± 0.11 mas 
along the major a xis, consistent w ith the size of the images of the first epoch as well as 
theoretical models (jHarmanedl2002l ) and the Ha disk of lSchmitt et al.l (120081 ) . The minor axis 



of th e disk is 0.63 ± 0.07 mas, larger than that expected in theoretical models (IBisikalo et al. 



2000l ; lLinnellll2000l ). implying this extended structure is perhaps from the electron scattering 



and /or free- free emissio n from the halo above the poles of the gainer (IJameson fc Longmore 



19761 : IZeilik et al.lll982l ). 



We also compared the flux ratios from our models with those obtained from H band 
light curves. The light curves from literatures^ give a value of = 1.86, where fx, /2 and 

f p are the fluxes of the donor, the disk of the gainer and the flux at the p rimary eclipse . 
Because the donor is not completely eclipsed by the disk at primary eclipse (ILinnelll 120001 ) . 
fp — f2 + a fii where a is the fraction of the donor flux that goes through. Therefore we can 
infer that jr > 0.86, consistent with our flux ratios derived from the models at phase 0.210 
and 0.828, i.e., 1.16 and 1.32. In addition, taking the average of the two values, ^ = 1.24, 
we get a = 0.165, suggesting that 16.5% of flux from the donor goes through at the primary 
eclipse. This also implies that the dono r contributes ~ 17% of total flux at the primary 
eclipses, consistent with the 20% value of ILinnelll (120001 ). 



The goodness of fit of the models (included in each panel of Figure [2]) are in general 



The light curve of lZeilik et all (| 19821) at the H band gives a —0.61 mag difference between the total 
m agnitude of the system and th e magnitude at primary eclipse. Interpolating the J and K band light curves 
of Jameson fc Longmore! (1976) we can obtain an H band magnitude difference of —0.74. Taking the average 
of the two we get a difference of 0.675 mag, corresponding to a ^^J 2 value of 1.86. 
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similar to that of the images. Nevertheless, discrepancies exist between the models and 
the images. For instance, the components' separations from the models are slightly smaller 
than the images. The gainer disk appears bigger in the models than in the images. These 
properties demonstrate the complexity of (3 Lyr and suggest that the disk may be asymmetric. 
They also highlight the deficiencies in our models and underscore the need for a self-consistent 
model that treats the radiative transfer properly and accounts for all epochs simultaneously. 



4. The orbit of (3 Lyr 



Despite the deficiencies in our simplified models, the positions of each component's 
center of light are well-constrained, especially near the maximum elongation of the orbit. 
Thus, the above results, along with the elements obtained from RV and light curve studies 
allow us to calculate the astrometric orbit of f3 Lyr for the first time. We ado pt P = 12. d 941 4 
and T mini = JD2454283.04 30 (on 2007Jul01) from the recent ephemeriai (Uk et alil^A 
and e = (jHarmaned 120021 ). The best-fit orbit using the model positions is shown in Figure 



|3j The resultant inclination, position angle of the ascending node (Q) and semi-major axis 
are listed in Table O Other orbital solutions using positions from the images are also listed 
in the table. Errors of orbital elements are estimated using Monte-Carlo simulations. The 3 
sets of inclination and Q in Table [3] are consistent with each other and suggest a retrograde 
orbit (i.e., position an gle decreasing with t ime). Our estimate on Q is roughly consistent 
with the 248?8 value in lSchmitt et al. ( 2008h , and is almost perpendicular to the orientation 
of the jet (163?5) implied bv lHoffman et al]Jl998h . 



We can also estimate the distance of (3 Lyr using orbital parallax (see Table [3]) by 
combining its angular semi-major axis with the linear asinz value, 57.87±0.62rJ| The 
distance from our models, 314±17 pc, is larger than that from the images, 278±24 pc and 
274 ±34 pc, but they a re all consistent within errors with the Hipparcos distance, 296±16 
pc ( Ivan Leeuwenl 120071 ). Finally, using the newly estimated Msin 3 i for both components 
(see footnote [3]) together with the inclination from the models, we get mass of the gainer = 
12.76±O.27M and mass of the donor= 2.83±O.18M . 



2 T m i„i (phase 0) is the epoch of primary minimum light. 

3 We obtain the semimajor axis by combining the semiamplitudc of the gainer K\ (41.4 ± 1.3 km s _1 , 
Harmanec & Scholz 1993; 42.1 ± 1.3 km s _1 [error assumed], Bisikalo et al. 2000; 35.4 ± 2.7 km s -1 [from 
their Fig. 5], Ak et al. 2007; yielding a weighted average of 41.1 ± 2.7 km s _1 ) with that for the donor K<z 
(185.27 ±0.20 km s -1 , Ak et al. 2007). We derive q = M 2 /M 1 = 0.222 ± 0.013, a sin i = 57.87 ± 0.62 R Q , 
Mi sin 3 i = 12.73 ± 0.27 M Q , and M 2 sin 3 i = 2.82 ± 0.18 M Q 
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5. Future work 

We have only presented simple two-component models in this work since we mostly 
focus on the orbital positions of ft Lyr. We have already discussed problematic discrepancies 
between the models and the images and also some internal inconsistencies between the model 
epochs. The systematic difference in component separations between the images compared to 
the model fits poses the most severe problem, limiting the accuracy of our distance estimates 
to ~ 15%. To address these issues and better understand other physical properties of ft 
Lyrae, a more physical, self-consistent model is required that treats the radiative transfer 
and the sizes of the two components properly, accounts for all epochs simultaneously, and 
incorporates the multi-wavelength information from eclipsing light curves. 
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04GI33G (JDM), and EU grant MOIF-CT-2004-002990 (NT). EP was formally supported by 
the Michelson Postdoctoral Fellowship and is currently supported by a Scottish Universities 
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Table 1. Observation logs for (3 Lyr 



Date (UT) 


Mean MJD 


Telescopes 






Calibrators 




2006Octl6 


54024.17 


W1-W2-S2-E2 


1 


29 Pe f 


■;, v And 




2007Jul03 


54284.25 


W1-W2-S1-E1 


3 


7 


Lyr, 


V 


Peg, v 


And 


2007Jul04 


54285.26 


W1-W2-S1-E1 


3 


7 


Lyr, 


V 


Peg, v 


And 


2007Jul07 


54288.22 


W1-W2-S1-E1 


3 


7 


Lyr, 


V 


Peg, a 


Cyg 


2007Jul09 


54290.25 


W1-W2-S1-E1 


3 


7 


Lyr, 


V 


Peg 




2007Jull2 


54293.26 


W1-W2-S1-E1 


3 


7 


Lyr, 


a 


Cyg 





a N(,(jfe : number of data blocks 

b Calibrator Diameters (milli-arcsec): 29 Peg = 1.017 ± 0.027, V And = 1.098 ± 
0.007 , a Cy g = 0.542 ± 0.021 (A . Merand 2008, private communication); 7 Lyr = 
0.74 ±0.10 jLegeett et al.lll986h ; v Peg = 1.01 ± 0.04 jBlackwell fc Lvnas-Gravl 
Il994h 



Table 2. Orbital positions of (3 Lyr 



MACIM BSMEM Model 



Date 


Phase 


Sep. 

(mas) 


P.A. 

(deg) 


Sep. 

(mas) 


P.A. 

(deg) 


Sep. 
(mas) 


P.A. 

(deg) 


Flux ratio 
(donor/gainer) 


2007Jul03 
2007Jul04 
2007Jul07 
2007Jul09 
2007Jull2 


0.132 
0.210 
0.438 
0.595 
0.828 


0.811 
0.891 

0.842 


255.4 
253.3 

72.3 


0.853 
0.886 

0.675 
0.783 


253.7 
254.4 

73.9 
69.6 


0.701 ± 0.091 
0.852 ± 0.045 
0.338 ± 0.105 
0.454 ± 0.042 
0.754 ± 0.063 


256.3 ± 4.0 
254.2 ± 2.1 
250.8 ± 7.3 
77.9 ± 1.4 
73.2 ± 0.8 


1.01 ± 0.11 

1 16 + - 20 
r- 10 -0.15 

3.51 ± 1.27 

2.43 ± 0.28 
i ,, +0.67 
-0.27 



a Some positions are omitted for images whose centroids cannot be separated. 
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Table 3. Parameters of (3 Lyr 





MACIM 


BSMEM 


Model 


Inclination (deg) 


92.10 ± 1.24 


91.96 ± 1.65 


92.25 ± 0.82 


n (dc g ) 


253.22 ± 1.97 


251.87 ± 1.83 


254.39 ± 0.83 


semi-major axis (mas) 


0.976 ± 0.083 


0.993 ± 0.122 


0.865 ± 0.048 


Distance (pc) 


276 ± 23 


271 ± 33 


312 ± 17 



400 



" li-wi 2007Jul04 



, S1-W2 
T E1-W1 
- □ E1-W2 

200 - G W1 " W2 



€ 

o 



-200- 



-400 L 



400 200 -200 -400 

East (m) 



Fig. 1. — Telescope spatial coverage of (3 Lyr on UT 2007Jul04, using the W1-W2-S1-E1 
configuration of CHARA. The symbols stand for different baselines. The longest projected 
baseline in this observation is 328.5m, corresponding to a resolution of 0.52 milli-arcseconds 
in the H band. The actual UV coverage is similar to this spatial coverage but each point 
spreads over 8 wavelength channels. 
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MACIM Image 



BSMEM Image 



2006Oct16 $ = 0,035 



0.5 - 
0.0 - 
-0.5 - 



X, = 0.67 



2007JUI03 = 0.132 




JG-1-06 



2007JUI04 $ = 0,210 




K = 0.42 



2007JUI07 <f = 0.438 




x; = 0.93 
2007JUI09 <f = 0.595 



0.5 - 
-0.0 - 
-0.5 




x;=o.8i 

2007JUI12 <f = 0.828 




K = 0.58 



X» = 0.28 




x: = o.3o 




x: = o.87 



x: = o.79 




X» = 0.52 



Model 




x;=i.42 {J. 



Too unresolved at 
primary eclipse for 
reliable modeling 




5C? = 131 



K = 0.46 



% = 0.69 



t = 0.66 



t = 0.52 



0.5 0.0 -0.5 -1 
East (mas) 



0.5 0.0 -0.5 -1.0 0.5 0.0 -0.5 
East (mas) East (mas) 



Fig. 2. — Reconstructed images and two- component models of (3 Lyr . The left, middle 
and right columns show the MACIM, BSMEM and model images respectively. Darker colors 
indicate higher intensity. The darker component is the donor. The contours in the images 
correspond to 0.3, 0.6, 0.9 of the peak intensity. Observing dates and corresponding phases 
(from the ephemeris in Ak et al. 2007) are labeled in the first column. The best-fit x 2 / DOF 
of each image is labeled in the bottom left corner. The resolution of the reconstructed images 
is 0.69 mas for the first epoch and 0.52 mas for the other 5 epochs, and the corresponding 
beams are shown in the first and second epochs in the middle panels respectively. Due to 
lack of enough resolution and the complexity of the radiative transfer at the first epoch when 
the star is behind the disk, no reliable model is available for our limited data. 
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0.4 


2007JUI12 




P Lyr orbit - 


0.2 








0.0 


7 2007JUI09 


• 




0.2 
0.4 






2007Jul07 ^^^^ ~ 
2007Jul04 ; 



1.0 0.5 0.0 -0.5 -1.0 



East (mas) 



Fig. 3. — The best fit relative orbit of (3 Lyr {solid line). The donor is indicated as a filled 
dot in the center. Positions of each epoch are shown by the open dots, surrounded by their 
error ellipses in dashed lines. The upper part of the ellipse is located towards the observer. 



